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Abstract 

Semileptonic and nonleptonic decays of the B c meson are studied in the frame- 
work of the relativistic quark model. The heavy quark expansion in inverse powers 
of the heavy (b, c) quark mass is used to simplify calculations while the light final 
quarks in the B and D mesons are treated relativistically. The decay form factors 
are explicitly expressed through the overlap integrals of the meson wave functions 
in the whole accessible kinematical range. The obtained results are compared with 
the predictions of other approaches. 

1 Introduction 

The investigation of weak decays of mesons composed of a heavy quark and antiquark gives 
a very important insight in the heavy quark dynamics. The properties of the B c meson 
are of special interest, since it is the only heavy meson consisting of two heavy quarks 
with different flavor. This difference of quark flavors forbids annihilation into gluons. As 
a result, the excited B c meson states lying below the BD production threshold undergo 
pionic or radiative transitions to the pseudoscalar ground state which is considerably more 
stable than corresponding charmonium or bottomonium states and decays only weakly. 
The Collider Detector at Fermilab (CDF) Collaboration [1] reported the discovery of the 
B c ground state in pp collisions. More experimental data are expected to come in near 
future from the Tevatron and Large Hadron Collider (LHC). 

The characteristic feature of the B c meson is that both quarks forming it are heavy and 
thus their weak decays give comparable contributions to the total decay rate. Therefore 
it is necessary to consider both the b quark decays b — > c,u with the c quark being 
a spectator and c quark decays c — *> s,d with b quark being a spectator. The former 
transitions lead to semileptonic decays to charmonium and D mesons while the latter 
lead to decays to B s and B mesons. The estimates of the B c decay rates indicate that 
the c quark decays give the dominant contribution (~ 70%) while the b quark decays 
and weak annihilation contribute about 20% and 10%, respectively (for a recent review 
see e.g. [2] and references therein). However, from the experimental point of view the 
B c decays to charmonium are easier to identify. Indeed, CDF observed B c mesons [1] 
analysing their semileptonic decays B c — > J/iplv. 

The important difference between the B c semileptonic decays induced by b — > c, u and 
c — > s,d transitions lies in the substantial difference of their kinematical ranges. In the 
case of B c decays to charmonium and mesons the kinematical range (the square of 
momentum transfer to the lepton pair varies from to g^ ax ~ 10 GeV 2 for decays to J /if} 
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and g max ~ 18 GeV 2 for decays to D mesons) is considerably broader than for decays to 
£?M an d _g(*) mesons (? max ~ 0.8 GeV 2 for decays to B s and g max 1 GeV 2 for decays 
to -B mesons). As a result in the _B C meson rest frame the maximum recoil momentum of 
the final charmonium and D mesons is of the same order of magnitude as their masses, 
while the maximum recoil momentum of the B^ and B^ mesons is considerably smaller 
than the meson masses. 

In this talk we consider semileptonic B c decays in the framework of the relativistic 
quark model based on the quasipotential approach in quantum field theory. This model 
has been successfully applied for the calculations of mass spectra, radiative and weak 
decays of heavy quarkonia and heavy-light mesons [3, 4, 5, 6, 7, 8]. In our recent paper 
[9] we applied this model for the investigation of properties of the B c meson and heavy 
quarkonia. The relativistic wave functions obtained there are used for the calculation 
of the transition matrix elements. The consistent theoretical description of B c decays 
requires a reliable determination of the q 2 dependence of the decay amplitudes in the 
whole kinematical range. In most previous calculations the corresponding decay form 
factors were determined only at one kinematical point either q 2 = or q 2 = q 2 nax and 
then extrapolated to the allowed kinematical range using some phenomenological ansatz 
(mainly (di)pole or Gaussian). Our aim is to explicitly determine the q 2 dependence of 
form factors in the whole kinematical range in order to avoid extrapolations thus reducing 
uncertainties. 



2 Relativistic quark model 

In the quasipotential approach a meson is described by the wave function of the bound 
quark- ant iquark state, which satisfies the quasipotential equation [10] of the Schrodinger 
type [11] 

(^-£)-*)-/^«io.«<„, d) 

where the relativistic reduced mass is 

E X E 2 M A -{m\-m 2 ) 2 

m ~ e^Te~ 2 - up ' (2) 

and Ei, E 2 are the center of mass energies on mass shell given by 

_ M 2 -m 2 + mj _ M 2 -m 2 + m 2 

1 2M ' 2 2M U 

Here M = E\ + E 2 is the meson mass, mi i2 are the quark masses, and p is their relative 
momentum. In the center of mass system the relative momentum squared on mass shell 
reads 

2 = [M 2 -(m 1 +m 2 ) 2 }[M 2 -(m 1 -m 2 ) 2 ] 
1 ] AM 2 ' 1 ' 

The kernel V(p, q; M) in Eq. (1) is the quasipotential operator of the quark- ant iquark 
interaction. It is constructed with the help of the off-mass-shell scattering amplitude, 
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projected onto the positive energy states. Constructing the quasipotential of the quark- 
antiquark interaction, we have assumed that the effective interaction is the sum of the 
usual one-gluon exchange term and the mixture of long-range vector and scalar linear 
confining potentials, where the vector confining potential contains the Pauli interaction. 
The quasipotential is then defined by [3] 

V(p, q; M) = «i(p)« 2 (-p)V(p, q; M) Ul (q)u 2 (-q), (5) 

with 

v( P , q; m) = ^zv(k) 7l ^ + Konf(k)r?r 2;/i + Konf(k), 

where a s is the QCD coupling constant, is the gluon propagator and k = p — q; 7 M 
and u{p) are the Dirac matrices and spinors. The effective long-range vector vertex is 
given by 

r„00 = 7m + ^<V*" » (6) 

where k is the Pauli interaction constant characterizing the long-range anomalous chromo- 
magnetic moment of quarks. Vector and scalar confining potentials in the nonrelativistic 
limit reduce to 

V v (r) = (1 - e)Ar + B, V s (r) = eAr, (7) 

reproducing 

Konf(r) = V s (r) + V v (r) = Ar + B, (8) 

where e is the mixing coefficient. 

The expression for the quasipotential of the heavy quarkonia, expanded in v 2 /c 2 with- 
out and with retardation corrections to the confining potential, can be found in Refs. [3] 
and [9, 4], respectively. The structure of the spin-dependent interaction is in agreement 
with the parameterization of Eichten and Feinberg [12]. The quasipotential for the in- 
teraction of a heavy quark with a light antiquark without employing the expansion in 
inverse powers of the light quark mass is given in Ref. [5]. All the parameters of our 
model like quark masses, parameters of the linear confining potential A and B, mixing 
coefficient e and anomalous chromomagnetic quark moment k are fixed from the analysis 
of heavy quarkonium masses [3] and radiative decays [6]. The quark masses m;, = 4.88 
GeV, m c = 1.55 GeV, m s = 0.50 GeV, m u ^ = 0.33 GeV and the parameters of the 
linear potential A = 0.18 GeV 2 and B = —0.16 GeV have usual values of quark models. 
The value of the mixing coefficient of vector and scalar confining potentials e — — 1 has 
been determined from the consideration of the heavy quark expansion for the semilep- 
tonic B — > D decays [7] and charmonium radiative decays [6]. Finally, the universal 
Pauli interaction constant k — —1 has been fixed from the analysis of the fine splitting 
of heavy quarkonia 3 Pj- states [3]. Note that the long-range magnetic contribution to 
the potential in our model is proportional to (1 + k) and thus vanishes for the chosen 
value of k — — 1. It has been known for a long time that the correct reproduction of the 
spin-dependent quark- antiquark interaction requires either assuming the scalar confine- 
ment or equivalently introducing the Pauli interaction with k — — 1 [13, 3, 4] in the vector 
confinement. 
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Figure 1: (a) Lowest order vertex function contributing to the current matrix element 
(9). (b) Vertex function taking the quark interaction into account. Dashed lines 
correspond to the effective potential V in (5). Bold lines denote the negative-energy part 
of the quark propagator. 

3 Matrix elements of the electroweak current 

In order to calculate the exclusive semileptonic decay rate of the B c meson, it is necessary 
to determine the corresponding matrix element of the weak current between meson states. 
In the quasipotential approach, the matrix element of the weak current J^ 7 = ?7 M (1— 75)^, 
associated with b — > q (Q = b and q = c, u) or c — * q (Q = c and q — s,d) transitions, 
between a B c meson with mass Mb c and momentum pb c and a final meson F (F = 
■0, D, B s , B) with mass M F and momentum p F takes the form [14] 

(F(p F )\J?\B c (p Bc )) = J ^^* FpF ( P )r M (p,q)^ cPflc (q), (9) 

where T At (p, q) is the two-particle vertex function and $m Pm are the meson (M = B c , F) 
wave functions projected onto the positive energy states of quarks and boosted to the 
moving reference frame with momentum p^. 

The contributions to T come from Figs, la, lb. The contribution r^ 2 ^ is the con- 
sequence of the projection onto the positive-energy states. Note that the form of the 
relativistic corrections resulting from the vertex function r( 2 ) is explicitly dependent on 
the Lorentz structure of the quark-antiquark interaction. In the leading order of the v 2 /c 2 
expansion for heavy quarkonia and in the heavy quark limit mj )C — * 00 for heavy-light 
mesons only contributes, while contributes already at the subleading order. The 
vertex functions look like 

r^(P,q) = u q (p q h^-l 5 )uQ(q Q )(2n) 3 5(p Qf - q Q 0, (10) 

and 

r<?>(p,q) = u q {p q )u Q >{PQ>){li»{l ~ 7i) 7-(^p^y 7i V(pQ^ - q Q 

+V(p Q < - q Q eg( ^ ( e fc j gQ) 7iV(l - lt)}u Q (q Q )u Q ,(q Q ,), (11) 
where the superscripts "(1)" and "(2)" correspond to Figs, la and lb, k = p q — A; k' = 
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q Q + A; A = Pf — p Bc ; 



\p) 



e(p) - (m 7 ° + 7°(7p)) 

Hp) 



Here [14] 



F i=l 



QQ,Q' 



ZQ,Q'{<l)j^±Y.^ ) tPBM\ 



and are three four-vectors given by 




M(E + M) 



pip 



I 



£ = + M 2 . 



The general structure of the current matrix element (9) is rather complicated, because 
it is necessary to integrate both with respect to d 3 p and d 3 q. The 5-function in the 
expression (10) for the vertex function permits to perform one of these integrations. 
As a result the contribution of to the current matrix element has the usual structure 
of an overlap integral of meson wave functions and can be calculated exactly (without 
employing any expansion) in the whole kinematical range, if the wave functions of the 
initial and final meson are known. The situation with the contribution r^ 2 ^ is different. 
Here, instead of a 5-function, we have a complicated structure, containing the potential 
of the gg-interaction in meson. Thus in the general case we cannot get rid of one of the 
integrations in the contribution of to the matrix element (9). Therefore, it is necessary 
to use some additional considerations in order to simplify calculations. The main idea is 
to expand the vertex function T^ 2 \ given by (11), in such a way that it will be possible 
to use the quasipotential equation (1) in order to perform one of the integrations in the 
current matrix element (9). 

B c — > r) c transitions. The natural expansion parameters for B c decays to char- 
monium are the active heavy b and c quark masses as well as the spectator c quark mass. 
We carry such an expansion up to the second order in the ratios of the relative quark 
momentum p and binding energy to the heavy quark masses mj, jC . It is important to take 
into account the fact that in the case of weak B c decays caused by b — > c, u quark transi- 
tion the kinematically allowed range is large (|A max | = (M| c — M F )/(2M Bc ) ~ 2.4 GeV 
for decays to charmonium and ~ 2.8 GeV for decays to D mesons). This means that the 
recoil momentum A of a final meson is large in comparison to the relative momentum 
p of quarks inside a meson (~ 0.5 GeV), being of the same order as the heavy quark 
mass almost in the whole kinematical range. Thus we do not use expansions in powers 
of \A\/m b:C or \A\/M F: but approximate in the expression (11) for T^ 2 ) the heavy quark 



energies e btC (p + A) = yml c + (p + A) 2 by e biC (A) = \Jml c + A 2 , which become inde- 

pendent of the quark relative momentum p. Making these replacements and expansions 
we see that it is possible to integrate the current matrix element (9) either with respect 
to d 3 p or d 3 q using the quasipotential equation (1). Performing integrations and taking 
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the sum of the contributions of and we get the expression for the current matrix 
element, which contains ordinary overlap integrals of meson wave functions and is valid 
in the whole kinematical range. Thus this matrix element can be easily calculated using 
numerical wave functions found in our meson mass spectrum analysis [9, 4]. 

B c — > transitions. In this case the heavy b quark undergoes the weak transition 
to the light u quark. The constituent u quark mass is of the same order of magnitude 
as the relative momentum and binding energy, thus we cannot apply the expansion in 
inverse powers of its mass. Nevertheless, taking into account the fact that the recoil 
momentum of the final meson in this decay is large almost in the whole kinematical range 
(as it was discussed above), we can neglect the relative momentum p of quarks inside a 
meson with respect to the large recoil momentum A. Thus in the region of large recoil 
(| A | ^> |p|) we can use the same expressions of the contribution to the current matrix 
element both for the B c — > and B c — > ip,r) c transitions. Moreover, the smallness of 
the r^ 2 ) contribution, which is proportional to the small binding energy, and its weak 
dependence on momentum transfer allows one to extrapolate these formulae to the whole 
kinematical range. As numerical estimates show, such extrapolation introduces only small 
uncertainties. 

B c — > B S (B)(*^ transitions. The heavy c quark undergoes the weak transition to 
the light s or d quark. The constituent s, d quark masses are of the same order of magni- 
tude as the relative momentum and binding energy, thus we cannot apply the expansion 
in inverse powers of their masses. The heavy quark expansion in l/m^ significantly 
simplifies the structure of the I^ 2 -* contribution to the decay matrix element, but the mo- 
mentum p dependence of the light quark energies e q {p) still prevents to perform one of the 
integrations. It is important to note that the kinematically allowed range for B c decays to 
B s and B meson is not large (|A max | = (M| c - M|)/(2M Sc ) ~ 0.8 GeV for decays to B s 
and ~ 0.9 GeV for decays to B mesons). This means that the recoil momentum A of a 
final meson is of the same order as the relative momentum p of quarks inside a heavy-light 
meson (~ 0.5 GeV) in the whole kinematical range. Taking also into account that the 
final B s and B mesons are weakly bound [5], we can replace the light quark energies by the 
center of mass energies on mass shell e q (p) — > E q — (Mp — ml + m 2 q )/(2M F ). We used such 
a substitution in our analysis of heavy-light meson mass spectra [5] which allowed us to 
treat the light quark relativistically without an unjustified expansion in inverse powers of 
its mass. Making these replacements and expansions we see that it is possible to integrate 
the current matrix element (9) either with respect to d 3 p or d 3 q using the quasipotential 
equation (1). Performing integrations and taking the sum of the contributions and 
we get the expression for the current matrix element, which contains ordinary overlap 
integrals of meson wave functions and is valid in the whole kinematical range. 

4 B c decay form factors and weak decay rates 

The matrix elements of the weak current J w for B c decays to pseudoscalar mesons P can 
be parametrized by two invariant form factors: 



(P(p F )\qYb\B c (p Bc )) = f + (q 2 ) p^+p F 




M 2 P 



+fo(q 2 ) 




M 2 P 



(12) 
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Table 1: Form factors of weak B c decays. 
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3.44 


0.49 
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0.35 
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q 2 = 


0.96 


0.80 


8.91 


0.72 


2.83 


1.06 


q 2 


= 


0.39 


0.39 


3.94 


0.42 


2.89 


0.20 



where q = pb c — Pf] Mb c is the B c meson mass and Mp is the pseudoscalar meson mass. 

The corresponding matrix elements for B c decays to vector mesons V are parametrized 
by four form factors 

(V(p F )\q^b\B(p Bc )) = M^+M v e ^ ae:PB ^ (13) 
(V{p F )Wl,b\B{p Bc )) = 2M v A (q 2 ) € ^q^+(M Bc + M v )A 1 (q 2 )(e^- e -^qA 



q- \ q< 



. . ■> < •«/ M 2 — M 2 

-Mq) 



• q 



M Bc + Mi 



v 



Pb c +Pf q2 1 



(14) 



where My and e M are the mass and polarization vector of the final vector meson. The 
following relations hold for the form factors at the maximum recoil point of the final 
meson (q 2 = 0) 

/+(0) = /o(0), 

Mr + M V Mr — Mv 

In the limit of vanishing lepton mass, the form factors fo and A do not contribute to 
the semileptonic decay rates. However, they contribute to nonleptonic decay rates in the 
factorization approximation. 

Now we apply the method for the calculation of the decay matrix elements described 
in the previous section. The explicit expressions for the form factors can be found in 
Refs. [17, 18]. Their calculated values are given in Table I. 
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Table 2: Semileptonic decay rates T (in 10 GeV) of B c decays. 

Decay our [19] [20] [21] [22] [23] [24] [25] [26] [27] 

B c ^r] c eu 5.9 14 11 11.1 14.2 2.1(6.9) 8.6 6.8 4.3 8.31 

B c ^rj' c eu 0.46 0.60 0.73 0.3 0.605 

B c -> J/ipev 17.7 33 28 30.2 34.4 21.6(48.3) 17.5 19.4 16.8 20.3 

B c ^i)'ev 0.44 1.94 1.45 1.7 0.186 

B c ^Dev 0.019 0.26 0.059 0.049 0.094 0.005(0.03) 0.001 0.0853 

5 C -> D*ev 0.11 0.49 0.27 0.192 0.269 0.12(0.5) 0.06 0.204 

B c ^B s ev 12 29 59 14.3 26.6 11.1(12.9) 15 12.3 11.75 26.8 

B c ^B*eu 25 37 65 50.4 44.0 33.5(37.0) 34 19.0 32.56 34.6 

B c ^Bev 0.6 2.1 4.9 1.14 2.30 0.9(1.0) 0.59 1.90 

B c -> B*ev 1.7 2.3 8.5 3.53 3.32 2.8(3.2) 2.44 2.34 



The differential semileptonic decay rates can be expressed in terms of the form factors 
as follows. 

(a) B c — > Peu decays 



G 2 F ^\V qQ \ 
24tt 3 



2\\2 



(b) B c — > Vez/ decays 



dT 

dg 2 



where Gp is the Fermi constant, V q Q is the CKM matrix element, 

A = lAI 



(Ml + Ml y - q 



2\2 



N 

The helicity amplitudes are given by 

2M Bc A 



Mil. + My 



V(q 2 ) T 



(M Bc + M v f 
2M Bc A 



2M v ^q I 



(M Bc + M v )(Ml -M 2 V - q^A^q 2 



AM 2 R A 2 A , 2 , 



(15) 



(16) 



(17) 



(18) 



M Bc + My 

The results of our calculations of the semileptonic and nonleptonic (in the factorization 
approximation) decay rates of B c meson are given in Tables 2, 3 in comparison with 
predictions of other approaches based on quark models [19, 21, 22, 24, 25, 27], QCD sum 
rules [20] and on the application of heavy quark symmetry relations [23, 26] to the quark 
model. Our predictions for the CKM favored semileptonic B c decays to charmonium 
ground states are almost 2 times smaller than those of QCD sum rules [20] and quark 
models [19, 21, 22], but agree with quark model results [24, 25, 26, 27]. Note that the 
ratios of the B c — > J/ipev to B c — > r\ c ev decay rates have close values in all approaches 
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Table 3: Nonleptonic B c decay rates T (in 10~ 15 GeV). 
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120a! 


242a! 


247a! 


141a| 


82.3a! 




143.6a! 


0.5a! 


5.3a! 


0.51a! 


1.65a! 


l.03a! 






0.7a! 


4.4a! 


1.40a! 


2.98a! 


l.28a! 






0.13a! 


4.8a! 


0.38a! 


1.45a! 


0.53a! 






3.4a! 


I3.lal 


4.50a! 


5.96a! 


4.56a! 







Bi 
Bt~ 
Bt- 

Bt- 
Bt- 
Bt~ 
Bt~ 
Bt 
Bt 
Bt 
Bt 

Bi- 

Bt~ 
Bt' 
Bt~ 
Bt- 
Bt 

Bi- 
Bt- 
Bi- 
Bt- 

Bt- 
Bt- 
Bt- 
Bt- 

Bt- 

Bt- 
Bt- 
Bt- 
Bt- 
Bt- 

Bt- 

Bt- 
Bt- 
Bt- 
Bt- 



VcP + 

> J/lpTT + 

J/ipK* + 

* ^^* + 
-> ^'^+ 

* ^'^*+ 

> B°K* + 

> B*°K + 

* B+K° 

> B+K*° 

> B*+R° 

^B+n° 
^B+p° 
-> B*+n° 



B* 



P 



ai ; 2 are the Wilson coefficients in the operator product expansion. 
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Table 4: Branching fractions (in %) of exclusive B c decays calculated for the fixed values 
of the B c lifetime r Bc = 0.46 ps, ai = 1.14 for 6-quark decays and ai = 1.20, a 2 = —0.317 
for c-quark decays. 



Decay 
6-quark decays 



Br 



B c - 


-> r] c eu 


0.42 


B c - 


-> rf c ev 


0.032 


B c - 


> J/ipev 


1.23 


B c - 


-> -^'ez/ 


0.031 


B c - 


-> £>ez/ 


0.013 


B c - 


-> D*eu 


0.037 



B- 
Bt 
Bt 

Bt 

Bt 

Bt 
Bt - 



Decay 
J/j/jK*- 



Br 

0.085 B+ 
0.21 B+ 
0.061 5+ 
0.16 £+ 

0.007 Bt 

0.011 5+ - 

0.005 Bt 

0.010 5+ - 



Decay 



Br 



> r)'ir + 


0.017 


> VcP + 


0.036 




0.011 




0.018 


i c K + 


0.001 




0.002 




0.001 




0.001 



B s ev 
B* s eu 



c-quark decays 
B c 
B c 

B c -> Bev 
B c -> 5*ez/ 

£+ -> 5 S 7T+ 

Bt - B sP + 

Bt -> b> + 



5 C + 
5 C + 



0.84 
1.75 
0.042 
0.12 
2.52 
1.41 

I. 61 

II. 1 
0.21 



B 



B 



> B S K* + 
* B*K + 
B+ -> £?°7r+ 
£+ -> B°p+ 
Bt -> 5*°7r+ 
5+ -> B*°p+ 



B°K* + 
B*°K+ 



0.003 S+ - 

0.11 S+ - 

0.10 Bt - 

0.13 5+ - 

0.026 Bt - 
0.68 Bt 
0.009 £+ 
0.004 
0.004 



c 

5+ 



c 



B *o^*+ oo33 

» fi+i? 0.24 

• B+K*° 0.09 

•B* + A' 0.11 

5*+^*° 0.84 

-> fi + 7T° 0.004 

-> fi+p° 0.005 

> B* + 7T° 0.001 

> 0.024 



except [23]. In the case of semileptonic decays to radially excited charmonium states our 
prediction for the decay to the pseudoscalar r)' c state is consistent with others, while the 
one for the decay to ip' is considerably smaller (with the exception of Ref. [27]). For the 
CKM suppressed semileptonic decays of B c to D mesons our results are in agreement with 
those of Ref. [23]. Our predictions for the CKM favored semileptonic B c decays to B^ 
are smaller than those of QCD sum rules [20] and quark models [19, 21, 22, 27], but agree 
with quark model results [23, 24, 25, 26]. For the CKM suppressed semileptonic decays 
of B c to fiW mesons our results are in agreement with the ones based on the application 
of heavy quark symmetry relations [23, 26] to the quark model. 

As one sees from Tables 2, 3 the theoretical predictions for B c weak decay rates 
differ substantially. Thus experimental measurements of corresponding decay rates can 
discriminate between various approaches. 



5 Conclusions 

In this talk we considered weak semileptonic and nonleptonic B c decays in the framework 
of the relativistic quark model based on the quasipotential approach in quantum field 
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theory. The weak decay form factors were calculated explicitly in the whole kinematical 
range using the heavy quark expansion for the initial active heavy quark Q = b,c and 
spectator antiquark Q' = c, b. The final light quark was treated completely relativistically 
without applying the unjustified expansion in inverse powers of its mass. The leading 
order contribution of the heavy quark expansion was treated exactly, while in calculating 
the subleading order contribution additional simplifying replacements were used. It was 
shown that such substitutions introduce only minor errors which are of the same order 
as the higher order terms in the heavy quark expansion. Thus the decay form factors 
were evaluated up to the subleading order of the heavy quark expansion. The overall 
subleading contributions are small and weakly depend on the momentum transfer q 2 . 

We calculated semileptonic and nonleptonic (in the factorization approximation) B c 
decay rates. Our predictions for the branching fractions are summarized in Table 4, 
where we use the central experimental value of the B c meson lifetime [16]. From this 
table we see that weak B c semileptonic decays to the ground and first radially excited 
states of charmonium and to D mesons, associated with b — > c, u quark transition, yield 
~ 1.7% and corresponding energetic nonleptonic decays (to charmonium and or n, p 
mesons) contribute ~ 0.6%. The semileptonic decays to B s and B mesons, associated with 
c — > s,d quark transition, give in total ~ 2.0% of the B c decay rate, while the energetic 
nonleptonic decays provide the dominant contribution ~ 19.3%. All these decays (to B s , 
B, charmonium and D mesons) add up to ~ 23.6% of the B c total decay rate. 

The authors express their gratitude to V. Kiselev, J. Korner, M. Muller-Preussker, V. 
Papadimitriou, V. Savrin and A. Vairo for support and discussions. Two of us (R.N.F and 
V.O.G.) were supported in part by the Deutsche Forschungsgemeinschaft under contract 
Eb 139/2-2. 
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